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Newly born neutron stars can present differential rotation, even if later it should be suppressed 
by viscosity or a sufficiently strong magnetic field. And in this early stage of its life, a neutron 
star is expected to have a strong emission of gravitational waves, which could be infiuenced by the 
differential rotation. We present here a new formalism for modelling differentially rotating neutron 
stars: working on the slow rotation approximation and assuming a small degree of differential 
rotation, we show that it is possible to separate variables in the Einstein field equations. The 
dragging of inertial frames is determined by solving three decoupled ODEs. After we establish our 
equilibrium model, we explore the infiuence of the differential rotation on the f and r-modes of 
oscillation of the neutron star in the Cowling approximation, and we also analyze an effect of the 
differential rotation on the emission of gravitational radiation from the f-modes. We see that the 
gravitational radiation from the f-modes is slightly suppressed by introducing differential rotation 
to the equilibrium stars. 



I. INTRODUCTION 

Differential rotation, until it becomes su_ppressed by 
viscosity or strong enough magnetic fields [J , might 
play an important role in the evolution of a newly born 
neutron star. (For a typical neutron star it takes between 
10-100 years to become uniformly rotating Q.) The equi- 
librium stellar models representing neutron star's differ- 
ential rotation were explored in some older papers 0, Q j 
and the oscillation frequencies for some types of fluid 
modes were calculated later in [1, . 

In this work we explore the evolution of linear pertur- 
bations in a slowly rotating neutron star with a poly- 
tropic equation of state. Moreover, the rotation profile 
of the star represents a first order deviation from the 
uniform rotation. We generalize the old semi-analytical 
results of Hartle for the uniformly rotating equilibrium 
model Q to a small deviation from the uniform rotation 
following the relativistic j-constant law. In particular we 
show that under a first order deviation from the slow uni- 
form rotation one can still separate spherical harmonics 
and obtain only a very small number of non-zero terms in 
their expansion. Furthermore, similar to the uniformly 
rotating case Q , one can find an exact analytic solution 
for the metric dragging function outside the star. 

In Section|lI]of this paper we use the consistency condi- 
tions imposed by the equilibrium model to constrain the 
value of the parameter representing the differential rota- 
tion, (for the given equilibrium parameters of the star). 
Such a relatively simple equilibrium model with the con- 
strained value of the differential rotation parameter is 
then used in Section lllli in the Cowling approximation. 
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to explore the various types of fluid modes. To obtain 
numerical results for the modes, in Section HVl we evolve 
the initial perturbation in time using a 2D Lax-Wendroff 
scheme [lO|. Our choice of a time evolution treatment 
is motivated by the co-rotation problem that affects the 
usual eigenvalue approach for obtaining the mode fre- 
quencies. The f-modes were computed for the slow rota- 
tion case in 0] and we compare our results with the 
results from those papers, further constraining the do- 
main of validity of the differential rotation parameter. In 
addition to the results of @, Q i we explore in more detail 
the behavior of the f-modes as a function of small values 
of the differential rotation parameter and we also find 
numerical values for some of the r-mode frequencies. (As 
can be seen in [ll[ the r-modes excite differential rotation 
in the stellar fluid.) Also, the f-mode eigenfunctions are 
extracted by using a pointwise discrete Fourier transform 
(DFT) on the evolution data. As the equilibrium stars 
have slow rotation with a low degree of differential rota- 
tion, the eigenmodes extracted show small change from 
their non- rotating counterpart. In Section |Vl by using a 
quadrupole estimate of gravitational emission timescale, 
we see a tendency that the differential rotation slightly 
suppresses the gravitational emission. These results are 
consistent with the results of [12] for the rapidly rotating 
stars. Finally, we finish with our conclusions in Section 
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II. EQUILLIBRIUM STELLAR MODEL 

Consider the background spacetime of a slowly rotat- 
ing star: 

ds^ = -e^dt^ -f e^dr^ -f r^dO^ ^ siTi^(e) [d</> - ujdtf , 

where u and A are functions of r and lo — uj(r,9) is the 
frame dragging function. We will use a polytropic equa- 
tion of state p = Xe^+^Z^, where p is the pressure and 
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e is the rest mass energy density of the star, 
rotation is described by the 4- velocity 



The fluid Further, we consider that the rotation of the fluid, f2, 
obeys the j-constant law: 



l + 7-r2 sin2(6')e-'' 



(1) 



The 7 parameter in the equation ([T]) describes the level 
of the differential rotation of the star. Then the TOV 
equations remain unchanged under the slow rotation (in 
the linear order in $7), and the frame dragging parameter 
uj has to be a solution of the equation: 
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l+7-r2sin2(e')e-'' 



= 0. (2) 



Furthermore, take the differential rotation to be also 
small, representing only a linear order perturbation from 
the uniform rotation case: < 7 << 1. Then one can 
expand a; in the 7 parameter as: 



uj(r, 0) = LUo{r, 0) + J ■ LOi{r, 0) + 0(72). 



(3) 
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Here corresponds to the case of uniform rotation and 
as we know from it depends only on r. (In the zeroth 
order expansion of ([2]) in 7 we obtain Hartle's equation 
for the uniform rotation problem Q.) Now, let us write 
the first order equation (in 7), which represents a cor- 
rection due to a small differential rotation modifying the 
uniform rotation problem. This will be: 



J__9 

^4 Qj. 



g(A-i/)/2 Q 

r2 s\n^{0)'d0 



sm^{0) 



duji 



167re(^-'')/2 . (e +p) [^^ _ (^^ _ ficy sin2(6')e-''] = 0. (4) 



Due to the fact that ujq does not depend on 0, one can 
simplify the problem by decomposing the terms in the 

I 



I 

equation ([4]) into the vector spherical harmonics and one 
obtains (symbol " ' " means r-derivative): 



r 4 • / 1 



j r^[{e{£ + 1) - 2}uu + W-rrr^e + p){uju - r'^e-'Ce[uQ ~ f^J}] 



(5) 



Here j{r) = e^^^^^^l"^ and Ct is a i — th coefficient of the 
decomposition of sin^{0) into vector spherical harmonics. 
One can express the decomposition as 



(6) 



and thus the only two non-zero coefficients Ci of the de- 
composintion are Ci =4/5 and C3 = —2/15. 

The two linearly independent solutions of the equation 
([S]) behave for £ > 1, both close to zero and at infinity, 
as 



(7) 



The £ = 1 case shows the same behavior ([7]) at the in- 
finity, but close to zero one has to be more careful: One 
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can try to Taylor expand the solutions at the origin, prov- 
ing that only one of the solutions is analytic around zero. 
Then one can naturally expect, that also in the case £ = 1 
one of the solutions is singular at zero. (For more details 
see Ell) 

It can be easily shown that in case Ci = 0, £ > 1, 
the regular behavior of the solution at the infinity can- 
not be matched with the regular behavior at zero (see 
again appendix lA II) . This means for Ci = no non- 
trivial relevant solutions exist. On the other hand, for 
Ci ^ we can, through the Green function, construct 
everywhere regular non-trivial solutions. This means the 
frame dragging function uj can be expressed in the linear 
order of 7 as: 

uj^ujo{r)+j-[ujn{r)+uji3{r){5cos^{0)-l}] . (8) 
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The fact that uJo{r) can be analytically solved outside 
the star is a known result and the solution is given 
as (r > R): 



Wo(r) 



(9) 



Moreover, similar to the uniformly rotating case, one can 
also find analytic solutions for wn, W13 outside the star. 
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(This is because the equation ([5]) can be rewritten outside 
the star into the form of the hypergeometric equation.) 
The analytic solutions are {r > R): 



and 



B2 



(10) 



Wl3(r) = B, 



3 • 



30 



210 - 180z + In 



z-1 



(120 - 300z + 180^2) 



(11) 



with z = r/2M. (M being the mass of the star.) Al- 
though this is maybe not obvious, the solution for W13 



can be shown to behave as 



when approaching in- 



finity. (All the terms with powers higher than cancel 
out. The derivation of the solutions and their asymptotic 
behavior is left for the appendix lA 21 ) Let us also add 
that the physical meaning of the constants _Bi,2 is the 
following: 



J, 



(12) 



where J is the angular momentum of the star. 

The procedure for numerically computing the frame 
dragging oj is the following: First we numerically deter- 
mine ojq by solving eq. (HJ inside the star for 7 = 6 (the 
usual Hartle equation from Q) with u}q{0) = and ujoiO) 
finite. The value iOo{R) then fixes the constant Bi from 
eq. ([9]), determining the behavior of ujq outside the star. 
The other components ujn and loi^ are obtained from the 
numerical integration of eq. (O with £ = 1 and 3, respec- 
tively. In order to pick the regular solutions, we write wn 
and Wis up to the second order in r with regular series 
expansions near the center as 



wii(r) = 60 



Btt 



and the constants bo, cq and B2 and B3 (see eqs. (fTOj) 
and (dJ)) are determined with a shooting method, by 
requiring that both cjn and ljis and their first derivatives 
be continuous on the stellar surface. Finally, the total 
dragging a;(r, 0) is computed with eq. ([5]). 

In the figures [T] and [21 we present some plots of the 
angular velocity f2 as a function of the radial coordinate, 
for different values of angle 6 and different values of the 
inverted 7 parameter. We also present (figures [3] and |31) 
plots of the frame dragging cj as a function of the radial 
coordinate, for different values of angle 6 and different 
values of the inverted 7 parameter. (The star is taken 
in the units c — G — Mq — 1 with the compactness 
M I R = 0.15 and with the equation of state parameters 
N = 1, K = IQO. This gives the stellar mass to be 
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M — 1.4. Note also that unless explicitly stated oth- 
erwise, we will use everywhere in the paper the units 
M© = 1.) In the figures [SJ [SI E] we show the 
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FIG. 1. Angular velocity Q, as a, function of r for different 
angles 9. (7 = 10"^) 

dependence of the central and the surface angular veloc- 
ities on the 7 parameter. (All the plots are taken at the 
equatorial plane.) The angular velocities are normalized 
by the Keplerian mass shedding limit, ^k- The minimal 
bounds on the 7 parameter are given by the equilibrium 
model, when either the central angular velocity reaches 
the value ^ O.SrJ/f , or when the surface angular velocity 
reaches zero. In this sense the minimal bounds on 7 are 
obtained naturally in the equatorial plane, as one can 
easily analytically observe that the second bound on 7, 
given by the surface angular velocity, has lowest value in 
the equatorial plane. (The first bound given by the cen- 
tral angular velocity is independent on 6.) In the next 
section, after computing the f-modes we further restrict 
the value of 7 by comparing our results for the f-mode 
frequencies with the results of Q- We confirm there is 
a very good agreement (less than 3 % error) up to the 
value 7 ~ 30~^, but for 7 ^ 20~^ the error is already 25 
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FIG. 2. Angular velocity f2 as a function of r for different 
values oi A = 7~^''^. 



fact that the perturbed fluid is barotropic: 
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FIG. 3. Frame dragging function as a function of r for differ- 
ent values of angle 6. (7 — 10^^.) 



FIG. 4. Frame dragging function as a function of r for differ- 
ent values oi A = 7"^'^^. 
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%, SO the bound (73) on 7 can be put as 7^ < 20 



III. LINEARIZED PERTURBATION 
EQUATIONS FOR THE FLUID 

We work in the Cowling approximation, thus we 
have only fluid perturbation variables, in particular: 
Se, dp, 6u^ . There are two more principles one uses 
to reduce the number of the variables to four: the four- 
velocity normalization condition 5(u^Ufj) — 0, and the 



FIG. 5. The solid lines represent the central angular velocities 
and the dashed lines the angular velocities at the surface as 
functions of 7, for different values of angular momenta. We 
use the star with the compactnes M /R = 0.1 and the equation 
of state with A'^ = 1.5, K = 10.86. (The total mass of the star 
is M = 1.47.) 

The remaining variables are: Su^,6u^,6u'^,SQ, with 
6Q = 5p/{p -f e). The dynamical equations are obtained 
from the three independent components of the perturbed 
Euler equation S{{Sj^ + u^u^)T.^'') — and the perturbed 
energy conservation equation 6{ukT!1'^) = 0. The final 
four equations for the linearized dynamics of the fluid 
can be written as: 
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FIG. 6. The solid lines represent the central angular veloc- 
ity and the dashed lines the angular velocity at the surface 
as functions of 7, for different stars with different angular 
momenta. The mass of the star decreases for different lines 
representing different cases from left to right. The compact- 
ness of the star is in all the three cases MjR = 0.1 and the 
parameters of the equation of state are (from right to left in 
the plot) iV = 0.5, 1, 1.5, K = 78106, 100, 10.86. 
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FIG. 7. The solid line represents the central angular veloc- 
ity and the dashed line the angular velocity at the surface as 
functions of 7, for different stars with different angular mo- 
menta. The compactness of the stars grows for different lines 
representing different cases from left to right and takes the 
values M/R = 0.1 ,0.15, 0.2. The star has the equation of 
state with A*' = 1 and K — 100. (The corresponding stellar 
masses are M = 1.06, 1.4, 1.62.) 
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5u% + n ■ 5u\ + 
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(15) 
(16) 



e + p 



e-^r^ &m^{e) • (17 - w) • (Ju^ + J< + 5u% + fe"*. 



e+p 



h' r + ^ r 2 

lL _|_ _ 

2 r 



^l.'5^^''-e'^^'-^-cot(0)-5w^ 
2 / e+p ^ ^ 



(17) 



IV. NUMERICAL RESULTS FOR THE MODES 

We used a 2D Lax-WandrofF scheme for solving the 
perturbation equations and the frequencies for the fluid 
modes were obtained through the Fourier transform of 
the time evolution of Sp at a given point inside the star. 
For the numerical integration we used the form of the 
equations ((T4)) . ((T5|) . ([16]), pTl) rewritten in the variables 
{5p, /*}, where /* is a momentum-like variable defined 
as /* = {p + e)Su'^ (similarly to what was done in [3l for 
Newtonian polytropes). We used symmetric boundary 
conditions for Sp (at the equatorial plane) for the f-modes 
and antisymetric boundary conditions for Sp for the r- 
modes. Also we used the regularity condition at both the 
radial center and the rotational axis. For the r-modes we 
used the initial data from and for the f-modes we 
used the initial value conditions from [l^ . 

In figure [H] we have a representative power spectrum 
obtained from our evolution data. One can see in this 
plot the correction to the rotational split added by the 
differential rotation. The relative heights of the peaks 
are rather arbitrary, and depend only on the initial data 
used. The width of the peaks is caused by the numerical 
dissipation of the algorithm used: combining that and 
our comparisons with values from the literature, we es- 
timate that our numerical error in the determination of 
the frequencies (see below) is within 3 %. 

We present in figure [H detailed results for the f-modes 
in the appropriate range of validity of 7. (As previously 
mentioned for 7 < 7b < 20~^ the agreement with the 
results of [zj is within less than 3 % error.) We present 
in this figure an equivalent to the correction to the fre- 
quency given by the rotational splitting of the f-modes in 
the uniformly rotating case, normalized now by the sur- 
face angular velocity of the star at the equatorial plane, 
fie , in units of fix ■ The two different data sequences cor- 
respond to constant J sequences of two polytropic stars 
with different polytropic indexes N, but approximately 
the same compactness. We can see that, for larger value 
of N, the correction starts with lower values, but grows 
faster with increasing differential rotation. This very fast 



growth shows a limitation of our first order treatment of 
the differential rotation. As seen in Q, when second or- 
der terms are takend into account, this growth becomes 
much less steep and much more " well-behaved" . 

We also present in tables U and [ll] results for the r- 
mode frequencies, which were not computed in We 
computed r-modes also for the uniformly rotating case, 
compared our results with |16l| a nd they agreed again 
with less than 3 % error. In [lg| they used decompo- 
sition of the perturbations into spherical harmonics and 
after truncating the coupled equations at £max, they time 
evolved the ID wave function to obtain the frequencies. 
Table U presents the values of the r-modes frequencies for 
a sequence of stars with constant angular velocity at the 
center flc and increasing differential rotation, while table 
nil presents the same results for a sequence of stars with 
constant angular momentum J. 

In the figure [TU] we present an equivalent to the cor- 
rection given in figure 1^1 but using now as central value 
the r-mode frequency for a uniformly rotating star, aro- 
We used the results from both tables HI and [Til in order to 
calculate these results. Note that the effect of the differ- 
ential rotation seems to be much weaker for the r-modes 
than for the f-modes (the scale of the vertical axis is now 
in Hz, and not in kHz as it was in figure [S]). 



A 


cr,.(Hz) 


n,(Hz) 


J 


1000 


491.55 


347 


0.2135 


500 


491.49 


347 


0.2134 


100 


490.78 


343 


0.2114 


50 


489.92 


331 


0.2050 


40 


489.57 


321 


0.2002 


30 


486.92 


302 


0.1898 


20 


446.90 


245 


0.1603 



TABLE I. The table of r-mode frequencies ar with constant 
angular velocity at the center (fic) and changing A = 7"^''^. 
(fie is surface angular velocity at the equatorial plane.) The 
star is taken with iV = 1, K = 100, M/R = 0.15, flc = 347 
Hz, M = 1.4. 
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A 


crr(Hz) 


J7c(Hz) 




2262 


491.55 


347 


347 


126.14 


492.47 


349 


346 


68.03 


495.19 


354 


346 


39.26 


511.02 


371 


343 


29.24 


549.59 


393 


339 


22.06 


589.43 


436 


331 


19.91 


602.76 


463 


327 



TABLE II. The table of r-mode frequencies ar with constant 
J and changing A = "/~^^^. The star is taken with TV = 
I, K = 100, M/R = 0.15, J = 0.2135, M = 1.4. 
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FIG. 8. The power spectrum obtained from our time evolution 
data for two stars with the same angular momentum J, with 
uniform (soHd line) and differential rotation with 7 = 0.003 
(dashed hue). Both stars have J — 0.2 and equation of state 
with N — 1.5, K — 10.86 (compactness M/R — 0.14 and mass 
M = 1.5). 



V. AN EFFECT OF DIFFERENTIAL 
ROTATION ON GRAVITATIONAL WAVE 
EMISSION FROM THE F-MODE 



FIG. 9. The rotational correction for the /--mode frequencies 
within the range of reliabihty of A = 7~^^^ for the stars with 
the equation of state with N=l, 1.5, K=100, 10.86, and with 
compactness M/R=0.15, 0.14, (and masses M= 1.4, 1.5) in 
sequences with constant J = 0.2. The quantity plotted is 
analogous to the quantities defined in [l^ for the uniformly 
rotating case, with a the frequency of the /_ for the given 
value of A, and ctq the frequency of /_ for the correspondent 
non-rotating star. 
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We here study an effect of differential rotation on grav- 
itational radiation by using a simple analysis. We use 
the Newtonian mass quadrupole formula to evaluate the 
gravitational wave emission. Luminosity of an eigen- 
mode is computed by using the eigenfunction extracted 
by DFT. The luminosity is the quadratic functional of 
the eigenfunction. On the other hand, we compute the 
kinetic energy of the eigenmode, which is also a quadratic 
functional of the eigenfunction. Taking the ratio of the 
luminosity and the energy, we obtain an inverse of the 
damping timescale of the eigenmode due to gravitational 
radiation. By comparing the timescale for different de- 
grees of differential rotation, we evaluate how differen- 
tial rotation affects gravitational emission from the eigen- 



FIG. 10. The equivalent of the rotational correction for the 
r-modes for sequences with constant J and sequences with 
constant angular velocity at the center Qc- This are the values 
from tables [H [Til so the star is taken with N — 1, K = 100, 
M/R = 0.15, (and M = 1.4). 



mode oscillation. 



1. Extracting eigenfunctions 

We extract the eigenfunction of the f-modes using the 
procedure described as follows. On each spatial grid 
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point, we performed a DFT of the physical variables to 
extract their power spectra. An eigenmode corresponds 
to a peak in the spectrum whose frequency is constant in 
space and is shared by different physical variables. We 
approximate the peak with a Lorentzian profile and ex- 
tract the central frequency, the peak amplitude, and its 
width. The collection of the amplitudes on each grid 
point gives the absolute value of the eigenmode excited 
in the simulation. We obtained sufficiently smooth eigen- 
function profiles for the f-modes, but we failed to extract 
higher order p-modes. It might be that we need to pre- 
pare initial data that contains the p-mode component 
with a larger amplitude than our current cases. In figure 
[TDwe present some typical results obtained for Sp, fr, fe 
and for the /+-mode which limits to the £ = m = 2 
f-mode in the non-rotating limit. 



6p f, 

HI H 




2e-08 4e-08 6e-0f 




FIG. 11. Eigenfunction of £ = |m| = 2 /+ mode for a differ- 
entially rotating star with A = 25.89. The upper half of the 
meridional section of the star is shown, with R the equatorial 
coordinate distance and z the coordinate distance parallel to 
the rotational axis. Up-Left: Sp. Up-right: fr. Down-left: 
fg. Down-right: The coordinate distance is normalized 
by the stellar radius. 



We also obtain the eigenfunctions for the /--mode, 
which shows qualitatively very similar functional profiles. 
Althogh we tried to extract p-mode eigenfunctions, but 
it was been hindered by numerical noise. 



2. Radiation timescale 

Formulae for the gravitational radiation from stellar 
oscillations are found in [l7] in the Newtonian limit. The 
energy of linear perturbation is defined as Eq.(15) there. 
Since our model is in the Cowling approximatiori, we 
have no gravitational perturbation (which is SU in '17'|). 
Therefore, we approximate the linear perturbation en- 
ergy as. 



J P Sv^SvldV, (18) 

where we assume the potential energy is regarded as the 
same size as the kinetic energy (i.e. we assume "equipar- 
tition" as in a simple oscillator). 

Energy loss by gravitational radiation is given by 



e>m 



-6DI 



dt ^ 



(19) 



where 



AttG {£+\)(e + 2) 
c2^+i<'(£-l)[(2£+l)!!]2^ 



(20) 



and 



J Sp r%* 



'dV. 



(21) 



Since the amplitude of the eigenmodes extracted from 
the linear evolution is arbitrary, we should study the 
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damping timescale Tg to characterize the efficiency of 
gravitational radiation, 



1 dE 
E dt 



(22) 



instead of an absolute amount of energy radiated. 

To compute Tg numerically, we need to define and com- 
pute the perturbed quantities appearing in the equations 
above. We need to have Sv"^ (perturbed 3- velocity), 5p 
(perturbed mass density) and dV (3- volume element). 
As for the density we used the perturbed rest mass 
density. The volume element is defined in the spatial 
hypersurface with t =const., where t is Schwarzschild 
time. The 3-metric 7^ is naturally chosen as 7 = 
diag(e^'*', r^, sin^ 0). Thus the corresponding volume 
element is dV — e^r^ sinddrdddip. 

I 



As for the 3-velocity perturbation, we adopted the def- 
inition below. The 3-velocity perturbation is expressed 
by the perturbed 4- velocity components as 



(u* 



(23) 



where u'^ is the 4- velocity. Then the components of the 
perturbed velocity (in coordinate basis) are expressed by 
our basic variables fi {i — r, 6, (p) as 



r 



(e -|-p)u* 



(e +p)u* 



(24) 



(25) 



= [1 + n{n - uj)r^ sin^ 



(e +p)u* ' 



(26) 



where /* (r, 9, ip) are defined as before as p = {e+p)Su^ 
We have used here 



and 



1 sin^ 9{n - w) 



(27) 



(28) 



Together with the equilibrium values of p, and metric 
coefficients, these perturbed variables are used to com- 
pute E, dE/dt and tqw 



Results 




0.004 



We compare the gravitational damping timescale tgw 
for different degrees of differential rotation parametrized 
by 7 for a fixed value of total angular momentum of the 
equilibrium star. In figure [T^] we plot sequences of tqw 
for the £ — \m\ = 2 f-modes with the angular momentum 
J = 0.2. 

In figure I12[ tqw is normalized by the corresponding 
timescale for a uniformly rotating star tgw,o- 7 = 
corresponds to the uniformly rotating model. A larger 
value of Tgw means smaller amount of gravitational ra- 
diation from the eigenmode. We see that the emissivity of 
gravitational radiation from each mode is reduced by in- 
troducing differential rotation. For the counter-rotating 
f-mode, this may be partly because the eigenfrequency 
is decreasing as we increase the degree of differential ro- 
tation. However it does not explain the increase of tgw 
for the prograde f-mode, whose frequency increases as we 



FIG. 12. Damping timescale of the eigenmode due to gravita- 
tional radiation, tqw, for £ = |m| — 2 f-modes. The sequences 
are obtained by fixing the angular momentum and increasing 
7, the degree of differential rotation. The timescale is normal- 
ized by that of the same eigenmode in a uniformly rotating 
star with the same angular momentum. The solid curve is for 
the /_ mode (counter-rotating mode) and the dashed one is 
for the /+ mode (prograde mode). 



increase the degree of differential rotation. For the /+ 
mode with differential rotation, the emissivity enhance- 
ment due to the increase of the frequency may be can- 
celed by a modification of the eigenfunction from that 
of the uniformly rotating case which reduces mass mul- 
tipoles. 
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VI. CONCLUSIONS 

In this work we dealt with a slowly rotating relativistic 
polytrope, such that has a nearly uniform rotation pro- 
file. We generalized the old result of , (for the equilib- 
rium model), for the first order deviations from the uni- 
form rotation. Similar to Q we are able to also provide 
an analytical solution for the metric in the exterior of the 
star. Furthermore, we used our equilibrium model to nu- 
merically compute (in the Cowling approximation) both 
f and r-mode frequencies. We also estimated the range 
of validity of our first order approach in the differential 
rotation parameter 7: We used the consistency condi- 
tions of the equilibrium model to constrain the domain 
of the 7 parameter, and some further restrictions were 
obtained by comparing our results for the f-modes with 
the known results in the literature. We provided detailed 
plots of the f-modes for different polytropes with different 
compactness/angular momenta and also provided some 
new results for the r-mode frequencies. By using a DFT 
we extracted the low order f-mode eigencfunctions from 
the evolution data. With the eignfrequencies and their 
eigenfunctions, the damping time of the oscillation due to 
gravitational radiation was estimated. Along the stellar 
models with a constant value of total angular momen- 
tum, we see a larger damping time as we increase the 7 
parameter to characterize the degree of differential rota- 
tion. This suggests that the inclusion of the differential 
rotation with our functional form tends to suppress the 
emission of gravitational wave for f-modes. 
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Appendix A: The equilibrium model 

1. Analysis of the regularity of solutions 

Let us show that if = 0, £ > 1, there does not 
exist an everywhere regular solution. In such case we are 
looking for a solution of the equation: 

[r'j ^'u]' = (Al) 
j [^(^+ 1) - 2 + 167rr2(e+p)] wi^. 

However, if a solution goes to zero at r — and in the 
same time at r 00 it must have at some point a lo- 
cal maximum/local minimum depending on whether it is 
approaching the zero at infinity from the negative values 
(close to the infinity the solution is negative), or from 
the positive values (close to the infinity the solution is 
positive). In case is close to infinity positive it must 
have a local maximum where is positive, in case is 



close to infinity negative, it must have somewhere a local 
minimum where ojif is negative. Take the equation (lAip 
at the point of such a local minimum/ maximum (call it 
Te). Then, due to the fact that the first derivative of 
vanishes at r^. the equation (|Aip can be written as: 

rb [^iy.e - (A2) 
j rl + 1) - 2 + 167rr2(e + p)] w^. 

But this equation cannot be fulfilled for a very simple 
reason: In case of a local maximum {aj"^}^^ < 0, and thus 
the left side of the equation must be negative, whereas 
the right side of the equation has always the same sign 
as a;i£, and lou is at the local maximum positive. (As we 
said before the local maximum is taken to be such that 
is at the local maximum positive.) In case of local 
minimum the opposite holds: {a;"^}^^ > and thus the 
left side of the equation must be positive, whereas the 
right side of the equation must be negative, because 
is at the local minimum negative. This argumentation 
means that ujxi for any £ ^ 1,3 must be represented by 
a trivial, zero solution. But what about wn and W13? 
(With wii regularity is not an issue, but one still requires 
that it shall have zero at infinity.) Here the situation 
is very different, since on the right side of the equation 
(|A2[) appears another term proportional to C3 and this 
prevents us from determining the sign of the right side of 
the equation at the extremum of wn, ^13. This means 
both wii, uj\z are in the game, and one has to proceed 
further in the analysis. 



Convergence of solutions near zero for 1=1 

Let us analyse the solutions in case £ = 1 near zero. 
For i — I the equation ([5]) turns close to to a more 
complicated equation: 

r + 4 - K.r = (A3) 

with K = 167r[e(0) -|-p(0)]. This equation seems not to 
have analytic solution, but one can verify that one and 
only one of the solutions is regular at by decomposing 
cjii into McLaurin series (or Taylor series at zero) 

00 

wii = ^ a„r" 

n=0 

and after doing this one obtains 

• fli = 0, 

• ari+2 = K ■ Qn/ [n^ -I- 5n -I- 4). 

This means there is only one solution that can be de- 
composed close to to McLaurin series and that strongly 
indicates that the other solution is singular at 0. Thus 
qualitatively the case ^ = 1 is the same than other, £ > 1 
cases. 
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2. The frame dragging function solutions outside 
the star 

Take the equation ^ outside the star. Consider that 
outside the star holds {M being the mass of the star): 




Then the equation ([5]) can be rewritten outside the star 
in the form 

r{r~2M)uj'lf+A{r-2M)uj[t,-[£{e+l)-2]uju = 0. (A5) 

After redefining the variable r = 2M ■ z and some alge- 
bras one can rewrite the equation (|A5[) in the following 
form: 

z{l ~ z)uu + 4(1 - z)uj'u + W + 1) - 'A^u = 0. (A6) 

By " ' " we mean here a z-derivative. Now consider that 
(IA6p is a hypergeometric equation with coefficients that 

I 

wi3 ^Di{\- "^z + zA +D2 



We can see that the first term is the divergent term. 
It is slightly less obvious that the other term is actually 
the convergent part of the solution behaving as z~^. One 
has to substitute z — e^^ and take the McLaurin series 
expansion of 

I 

D2 ■ (-6^ - Se^ - 30e + 210 - ISOe"^ + 



J 



Wl3 = D2 
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I 

Then if one expresses the second solution term in e: 



• {l20-300e"^ -h 180e"^}), (AlO) 

I 
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terms up to the 5-th power cancel. This means ^13 term 
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